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1. Introduction 

The search for physics beyond the Standard Model is one of the main objectives of 
the upcoming experiments at the CERN Large Hadron Collider (LHC) as well as 
the currently running ones at the Fermilab Tevatron. Of special interest is the scalar 
(Higgs) sector, which presumably is responsible for electroweak symmetry breaking 
and generation of particle masses. The discovery of a neutral scalar Higgs boson (h°) 
would be a big step in understanding electroweak symmetry breaking. But at the 
same time it may be difficult to decide whether it is the Standard Model Higgs boson 
or if it belongs to for example a supersymmetric theory. However, the discovery of 
a charged Higgs boson would be a very clear signal of physics beyond the Standard 
Model, and would give valuable insight into the parameter space of this physics, 
whether supersymmetry or some other two Higgs doublet model. 

Present model-independent limits on the charged Higgs boson mass from LEP 
are > 78.6 GeV (95% CL) |T| whereas the limits from the Fermilab Tevatron 
are close to m H ± > 160 (130) GeV g | for large (small) tan/5 (tan/3 > 100 and 
tan/3 < 1 respectively). As usual, tan/5 = — is the ratio of the vacuum expectation 

values of the two Higgs doublets which determines the coupling strength to the top 
and bottom quarks (roughly Xn±tb ~ m b tan (3 + m t cot/3). For a recent review of 
the prospects of further charged Higgs boson searches at the Tevatron and the LHC 
see §. 
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Figure 1: Illustration of the relation between the gb — > H^t and 55 — > H^tb processes. 



At hadron colliders, the main contribution to direct single charged Higgs boson 
production is through the two twin-processes gb — > H^t and gg — > H^tb 1 . The 
reason for calling them twin-processes is simply that they describe the same under- 
lying physical process, as illustrated in fig. |], but using different approximations. 
In fact, the latter process enters in the next-to-leading order (NLO) calculation of 
the first one |5j, [/J. However, in the following we will concentrate on the simu- 
lation of charged Higgs boson production in Monte Carlo event generators such as 
Pythia m and Herwig || where only tree-processes are included 2 . Note that al- 
though our matching procedure works for all charged Higgs masses, we are mainly 
interested in the region m#± > m t —mb where the cross-sections for the two processes 
are of similar size. 

In the gb -> HH process the 6-quark is considered as a parton in the proton 
with a corresponding parton density which resums all the, potentially large, leading 
logs of the type (a s log(//F/raf,)) n that arise when integrating the DGLAP evolution 
equation from the threshold given by the 6-quark mass to the factorization scale //p, 
where the proton is probed. Naturally, the treatment of the 6-quark as a collinear 
massless parton in the proton leads to certain approximations which are at best only 
valid in the collinear limit. In contrast, the 6-quark is not considered as a massless 
parton in the gg — > H ± tb process and consequently the kinematics of the process are 
exact. This also means that the finite parts, which are not logarithmically enhanced, 
are also included correctly to order a 2 s . On the other hand the gg — > H ± tb process 
does not contain the (a s log(/ip/mb)) n terms with n > 1 which are resummed in the 
6-quark density. 

In case of the total cross-section it is straightforward how to combine the gb — > 
HH and gg -> HHb processes without introducing double counting |TTJ (see also |L2], 



T3j). One simply adds the cross- sect ions from the two processes and subtracts 



the common term, which is given by the gb' -> HH process where only the lead- 



^or brevity we make no distinction between quarks and anti-quarks unless it is not clear from 
charge-conservation what the correct assignments are. 

2 In the last years there has also been a development of new techniques for including NLO matrix 



elements in Monte Carlo generators, known as "MC@NLO" 1 1 . We will comment more on the 
relation of our method to MC@NLO below. 
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ing order logarithmically enhanced contribution to the 6-quark density, b'(x,{ip) oc 
a s \og(fj,F / nib) is used. In this paper we will show how this subtraction procedure 
can be extended also to the differential cross-section (irrespective of whether the 
outgoing 6-quark is observed or not). One has to keep in mind that even though 
the 6-quark is considered to be collinear in the gb — > H t process and thereby the 
accompanying 6-quark from the gluon splitting is also collinear, this is no longer true 
in a Monte Carlo event generator. The reason is that in the Monte Carlo, the hard 
process is combined with an initial state parton shower which "undoes" the DGLAP 
evolution back to a starting scale Q < m b , and thereby generates an accompanying 
6-quark with non-zero p± (transverse momentum with respect to the beam axis). 

The experimental techniques for detecting charged Higgs bosons rely heavily 
on tagging of 6-jets and/or hadronic r-jets together with missing p±. The 6-jets 
can originate from the charged Higgs boson decay, the top quark decay or from 
an accompanying 6-quark whereas the r-jets and missing p±_ of interest come from 
the charged Higgs boson decay. Below the top quark mass the dominating charged 
Higgs boson decay is H ± — ■> ru T whereas for heavier masses the H — > tb decay mode 
dominates. Thus the minimum requirement is to tag one 6-jet and either a r-jet 
together with missing p± or two additional 6-jets depending on the decay mode of 
the charged Higgs boson. In addition the accompanying 6-quark may also be tagged. 



Typical cuts used in studies by the ATLAS [14], [L5| and CMS JT6| collaborations for 
6-jets, r-jets and missing pj_ are given by the following: p± t b > 30 GeV, |^| < 2.5, 
P±,t > 100 GeV, \t] t \ < 2.5, and p_L im i SS > 100 GeV. With these cuts the 6-tagging 
efficiency is of the order of 50% with a miss-tagging rate of about 1% whereas the r- 
tagging efficiency is of the order of 40% jlfj. In addition the difference in polarization 
between a r-lepton coming from the scalar Higgs decay and the vector H^-boson decay 
leads to a harder spectrum when the r decays into one charged pion which in turn 
can be used to enhance the signal |17| . 

Conventionally the strategy, when investigating the prospects of detecting charged 
Higgs bosons at the LHC, has been to use the gb — > H ± t process combined with par- 
ton showers when the accompanying 6-quark is not tagged in the final state (see 
e.g. [|14], |U|) whereas the gg — ■> H^tb process has been used when the 6-quark is 
tagged (see e.g. [T^| and |T^]). As we will show in this paper, the only way to get 
reliable predictions in the latter case is to make a proper matching of the two pro- 
cesses, whereas in the first case the need for matching depends on the details of the 
selection used in defining the signal and the precision that one is aiming for. We 
will also see that in general it is not possible to simply divide the phase-space into 
two different parts, based for example on the p± of the accompanying 6-quark (as in 
e.g. pOl), with each of them only populated by one of the processes. 

The method that we present is general. The actual implementation has been 
done in Pythia although the same procedure could also be applied to Herwig. 

The outline of the rest of the paper is as follows. We start by recalling the proper 
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matching procedure in case of the total cross-section. In section 3 we then show how 
to generalize the method to be applicable also for differential cross-sections. The 
results of the matching procedure are presented in section 4 and in section 5 we 
discuss how to choose a proper factorization scale. Finally, in section 6 we give our 
conclusions and a short outlook. 



2. The total cross-section 



As already discussed in the introduction, the two different approximations describing 
H ± production at hadron colliders that we are interested in are, 

gb(b) -> tH~ (iH + ) (2.1) 
gg -> tbH~ (ibH + ) . (2.2) 

The first one (|2.1|) , which we will denote the leading order (LO) or 2 — ► 2 process, 
includes the logarithmic DGLAP resummation of gluon splitting to bb pairs via the 
6-quark density, b(/Mp) ~ ^](a s log(/ii?/m{,)) n , whereas the second one, which we will 
denote the 2 —* 3 process, retains the correct treatment of the accompanying 6-quark 
to order a 2 . 

In case of the total cross-section, the two approximations can be combined by 
simply adding them and subtracting the common term. The double-counting term 
which needs to be subtracted is, to leading logarithmic accuracy, given by jll 



o"dc = / dx 1 dx 2 



g{x u /j, 2 F )b'(x 2 , f*F) fa 2 fa (zx,x 2 ) + x 1 <-+ x 2 



(2.3) 



where b'(x,fi 2 F ) is the leading logarithmic contribution to the 6-quark density, 

A) = ^ lo S % J f W*) 9 (? ^) (2 ' 4) 

and Pg^ q q(z) = \ [z 2 + (1 — z) 2 ] is the splitting function. 
The matched integrated cross-section is thus given by 

a = a 2 ^ 2 + o"2^3 - odc- (2.5) 

The resulting cross-section for the process pp —> H^tb at the LHC is illustrated in 
fig. § as a function of m H ± for the case tan/3 = 30 and using fip = {m t + m H ±)/A. 
(The choice of factorization scale will be discussed in section |.) We have used a 
running 6-quark mass in the Yukawa coupling evaluated at , the renormalization 
scale for a s has also been set to m#± and we have used the CTEQ5L |21] parton 



densities. As can be seen from fig. |2], for charged Higgs masses below m t — rrib 
the 2^3 process dominates the cross-section. In this region, the 2^3 process 
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Figure 2: Integrated cross-section (components and matched total, using the exact 
double-counting term) as a function of the H mass at LHC. Here tan/? = 30 and \xp = 
(mt +m#±)/4. Note that the double-counting term contribution (DC) is subtracted from 
the sum and that the two different versions of the double-counting term corresponds to 
whether the kinematic constraints are included or not. 



can be well approximated by intermediate top production, gg — > ti — > tbH ± . This 
approximation breaks down for m H ± close to and above the top mass, see e.g. P^ ]. 
The matching procedure described in this paper, although it works for all charged 
Higgs masses, is of greatest interest in the regions where the 2^2 and 2 — > 3 
processes are of similar size, i.e. for m#± > m t — m&. At the same time our procedure 
gives a smooth transition between the light and heavy Higgs mass regions, which may 
be helpful when devising search strategies in this so called transition region. 

The way of defining the total cross-section in eq. ( |2.5| ) is implicitly assuming that 
the accompanying 6-quark from gluon splitting is observed (which at least in principle 
is always possible since there are no 6-quarks in the initial state). Thus we are in 
effect talking about a leading order (in a s ) cross-section for the process pp — > H ± tb. 
Using the power counting rules of ref. ||23|| , the first term in eq. (|2.5|) , o"2^2, is of order 
a 2 s log(/ii?/m(,) whereas the correction (02^3 — odc) is of order a 2 s . This 1/ log(/ii?/mh) 
correction, which arises from gluon splitting into a 66-pair, is also part of the NLO 
calculation of the cross-section for the pp — > H^t process. In addition the NLO 
calculation also includes virtual corrections as well as real corrections from gluon 
emission both of which are of order log(/ii?/mb). Another essential difference is of 
course that in the NLO calculation the accompanying 6-quark is not observable. 

Comparing with the MC@NLO approach as it is has been implemented for heavy 
flavor production P^| , this uses the equivalent of o 2^3 together with its NLO a s 
corrections as starting point. In other words, MC@NLO has so far only been based 
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on the so called flavor creation processes where there are no heavy quarks in the initial 
state and thus the equivalent of the 0^2 process is not included. As a consequence 
the MC@NLO method has not yet been applied to charged Higgs bosons production. 

A more complete expression for b', which also takes into account non-logarithmic 
contributions arising from kinematic constraints due to the finite center of mass 
energy and the non-zero 6-quark mass, will be derived in the next section. It is 
important to include these corrections also in the double-counting term since they 
are included both in the parton shower and the 2^3 matrix element. The difference 
when including the finite corrections are substantial, leading to a large reduction of 
the double-counting term, as can be seen from fig. |^. For example, in the case 
m H ± = 250 GeV, tan/3 = 30, and fip — (jn t + m H ±)/A the double-counting term is 
reduced from 0.81 to 0.55 pb. 

For the purpose of the generalization of this method to the differential cross- 
section, which will be given in the next section, it is instructive to reconsider the 
derivation of the integrated double-counting term. The starting point is the DGLAP 
evolution equation for a massless 6-quark which we write, 



Since we are only interested in the leading logarithmic contribution to the 6-quark 
density and it is generally assumed that the 6-quark density is only perturbatively 
generated 3 the last term can be dropped. Furthermore, in the massive case, the 1/k 2 
term in ( |2.6| ), which originates from the fe-quark propagator (the quark line marked 
with a dash in fig. [I]), has to be replaced with, 



1 1 1 



k 2 k 2 — ml Q 2 + ml 
where Q 2 = —k 2 . In summary we get the following expression, 



J2.7) 



where the last step is valid to leading order in a s . Eq. (|2.4j) is then obtained by 
simply integrating between the integration limits Qmin = m \ an d Qmax = V% which 
are correct to leading logarithmic accuracy. In the next section we will derive the 
exact integration limits based on kinematic constraints that have to be taken into 
account in the differential cross-section. 



3 Thus we assume that the intrinsic (2f| contribution to the 6-quark density can be neglected for 
our purposes. 
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3. Matching the differential cross-section 

In this section we give the details of how to extend the matching procedure for the 
total cross-section outlined above to also be valid for the differential cross-section. 
When doing this it is important to keep in mind that in a Monte Carlo generator 
such as Pythia, the leading order 2^2 process is supplemented with initial state 
parton showers which generate a non-zero p±^ for the outgoing 6-quark from the 
gluon splitting, based on the leading order DGLAP evolution equations. The p± } b~ 
distribution generated in this way is essentially of the type da/dp± ! b oc pi^/m^, 
where m\ b = ml + p\ b is the transverse mass squared, up to the maximal p± gen- 
erated by the parton shower, which in general 4 should be given by the factorization 
scale fip. This means that from the Monte Carlo point of view the final state parti- 
cles of the two processes (|2.1| ) and ( |2.2|) are the same, but generated with different 
distributions due to the different approximations used. 

In order to distinguish the two processes and define the double-counting term 
we need some "guiding principle" which tells us in which region of phase space the 
different approximations are valid. In our case it is natural to use the ^-distribution 
of the outgoing 6-quark. To illustrate the logic behind this choice it is instructive to 
consider the pj^-distribution from the matrix element describing the 2 — > 3 process, 
which is shown in fig. []. In order to get a better illustration of the behavior at small 
P_L,b the differential cross-section has been multiplied by m 2 _ L b /p± t b- 

From the figure it is clear that for small p±^ the differential cross-section for the 
2 — > 3 matrix element behaves more or less as da/dp± t b oc pj^/m^ b (similarly to the 
parton shower), whereas for large p±$ (where px,b is the hard scale of the process) it 
behaves as da/dp^b oc l/p 3 ±b , as seen by simple dimensional analysis. This means 
that, as one expects on general grounds, for small p± t b the 2^3 matrix element can 
be factorized into two parts, the g — > bb splitting and the gb — > tH^ hard process 
(or symbolically 2^3 = 1 ^ 2 ® 2 ^ 2). For large p±^ on the other hand, where 
the outgoing 6-quark is part of the hard process, it is not possible to factorize the 
process and one has to retain the complete matrix element. 

These observations can now be used to identify the regions of phase space where 
the two approximations have their respective strength. For small p±^ we are in 
the collinear region where the cross-section is dominated by logarithms of the type 
{a s \og(m\ b /m b )y i /i r n 2 Lb (the 2^3 matrix element only contains the leading 
term n = 1) and therefore resummation effects, which are included in the DGLAP 
evolution of the 6-quark density, are important. For large p± t b on the other hand we 
are in the hard region where factorization breaks down and we have to use the 2 — > 3 
matrix element. 

4 In Pythia the maximal p± in the parton shower is by default equal to the factorization scale 
Hp but it can also be set to be some factor times the factorization scale jip- 
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Figure 3: The differential cross-section da/dp± t b multiplied by (m 2 + p 2 L b)/P±,b f° r ^ e 
2^3 matrix element and the differential double-counting term, with [ip = (mnt + m _ff±)/4- 
Note the two distinct regions in the 2^3 cross-section, for small and large p± t b respectively. 



In order to obtain the correct differential cross-sections for the accompanying 
6-quark, we need to employ a matching procedure between the gb — > H^t and gg — > 
H^tb processes. Just as in the case of the total cross-section the matching is done 
by adding the two processes and subtracting the double-counting thus introduced in 
the collinear region. The most basic requirements of such a procedure are that all 
resulting differential cross-sections must be smooth, and that the integrated cross- 
section must be the correct one as given by eq. ( |2.5|) . Also the resulting differential 
cross-sections should be given by the leading order process for small p± t b and by the 
2 — > 3 matrix element for large p± t b- 

The key observation in the generalization of the matching procedure to the dif- 
ferential cross-section is that the leading logarithmic 6-quark density (|2.9| ) defines a 
differential distribution in the variables z and Q 2 , which together with the differen- 
tial cross-section for the 2 — ^ 2 process and two azimuthal angles gives a complete 
differential definition of the double-counting term. By picking events from this distri- 
bution and subtracting them in the data analysis (e.g. histograms), we can account 
for the double-counting not only for the total cross-section but also differentially. 5 

In this way, we can write down the complete expression for the double-counting 
term based on eqs. ( |2.3| ) and Q2.9| ) (properly accounting for the phase-space limita- 
tions as well as the mass of the incoming 6-quark): 



5 The double-counting term is by construction always smaller than the 2 — > 2 term in the whole 
phase-space, since b'(x,Q 2 ) < b(x,Q 2 ). Therefore, if only the number of events is large enough, 
there should be no risk that the matched cross-section turns out to be negative in any phase-space 
bin. 
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a DC = / — / dy*- ^d{cos9)\M h0 \ -^f^x 



2 ' 



T Ji\ OKT S 2 27T 

(3.1) 



) 2 

& max 



<i(<5 2 ) -''i .-n 2 „ , 2 



dzP g ^ qq (z) / — 2 — g(— , n F ) x 2 g(x 2 , n F ) + x x «-> x 2 

J QLn ^ + m b z z 



Here .Mlo is the matrix element for the leading order process ( |2.1|) and the variables 
for the 2 — > 2 process are defined as follows: r = xix 2 , x lj2 = y/re ±y * , s = ts, 
9 is the polar angle of the t-quark in the CM system of the 2 — > 2 scattering, and 
#34 = s -: V(s - m'i - m 2 H± ) 2 - 4m 2 m 2 H± . 
The integration limits are given by 

(3.2a) 



Train 


im t + m H ±) 2 /s 


^max 




(Q 2 op t + s)(Q 2 opt + m 2 ) 


O 2 ■ = 

°0 min 


\[s{z^-l)-ml\- 


c 2 = 

V max 


min \ [s (z _1 - 1 



> <5o P t = mm ( \/ smg, /4 ) (3.2b) 



s {z- 1 - 1) - m 2 ] 2 - Asm 2 (3.2c) 



wis] + iv[ s " x) " " 4Sm 'j • (3 - 2d) 

In deriving these limits we have identified z as the ratio of the center of mass energies 
of the 2^2 and 2^3 processes, z = s/s'. Assuming that the outgoing b- and 
i-quarks as well as the charged Higgs boson are on-shell, this identification of z gives 
the following expression for the transverse momentum of the outgoing 6-quark in the 
center of mass system of the two gluons: 

p2 2 _ 2 (£±0!)M±«!). ,3.3) 

s 

Note that Q 2 = —k 2 , where k is the 4-momentum of the 6-quark propagator. From 
this the limits on z and Q 2 of eq. (|3.2| ) follow by considering the collinear situation 
when p\ b = and taking into account that the virtuality cannot be larger than the 
factorization scale fi 2 F . In case of the upper limit on z one also has to find the optimal 
Q 2 value which maximizes z by setting dz/dQ 2 = 0. 

It should be noted that the renormalization scale (1r in eq. fl3.1| ) is the same 



as in the hard 2^2 subprocess (|2.1|) and that the gluon density related to the 
g — > bb splitting is evaluated at the same factorization scale fip as the fr'-density and 
the other incoming gluon. These choices corresponds to requiring that the double- 
counting term approximates the 2 — > 3 matrix element for small p± ) b- In the Pythia 
parton shower the renormalization scale is given by (1 — z)Q 2 ~ p 2 ± b and the factor- 
ization scale by Q 2 , and since the parton shower corresponds to a leading- logarithmic 
resummation, a s is calculated at 1-loop. For pj_b ~ /Af, the difference between the 
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different scale choices is negligible and thus the double-counting term will approach 
the distribution obtained from the 2 — > 2 process with parton showers in this region, 
ensuring that the matched differential cross-section smoothly interpolates between 
the 2 — > 2 process for small p±^ and the 2 — > 3 process for large p± t b- 
Taking the limit s — > oo, eq. ( |3.3| ) also gives the relation 

= ^ (3-4) 



Q 2 + m\ p\ h + m b 

which explains the plateau for small p±^ seen in fig. ^ for both the double-counting 
term and the 2^3 matrix element. The figure also illustrates the importance of 
the kinematic constraints in the double-counting term which otherwise would have 
been flat all the way up to p^,b = fJ-F- In section || we will discuss how the extension 
of this plateau can be used to set a proper factorization scale for the 2 — > 2 process. 

The differential cross-section for the double- counting term has been implemented 
as an external process to Pythia and will be made available for download p6| 
whereas the 2^3 process is implemented ctS ct TQ gular process in Pythia from 
version 6.223. 



4. Results of the matching procedure 

In this section we present the results from a case study of our matching procedure 
using m H ± = 250 GeV and tan/3 = 30 at LHC energies. 

We use Pythia, including initial- and final-state radiation, with default settings, 
except that the factorization scale is set to be /i^ = (m t + rriH±) / '4 (this choice of 
factorization scale will be motivated in section ||), and a s in the hard process and the 
running 6-quark mass entering the Yukawa coupling are evaluated at m#± (we use 
a 1-loop expression for a s and the Yukawa coupling since we are effectively dealing 
with a leading order process as argued in connection with eq. |2.5| ). In addition the 
outgoing 6-quark in the initial state parton shower as well as in the double-counting 
term is put on-shell in order to facilitate a more clearcut comparison with the 2^3 
process, where the outgoing 6-quark is on-shell. In a realistic case, using jet-finding 
algorithms, there should be no observable difference compared to having the outgoing 
6-quark off-shell in the parton shower. Finally the maximal scale for the initial state 
shower is set to the factorization scale and we use the CTEQ5L parton densities. 

Since we are only dealing with leading order calculations the overall normaliza- 
tion is quite uncertain, especially given that the cross-section is proportional to the 
renormalized 6-quark mass from the Yukawa coupling. Therefore, in the following 
we will concentrate on the shapes of the distributions. The results for differential 
cross-sections in a few different variables are presented in figs. [| and |5|. 

First of all it is clear from fig. ^a that the matched ^^-distribution looks as 
expected. It follows the leading order cross-section for small p± t b, where the 6-quark 
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Figure 4: Differential distributions in (a) p±^ and (b) r]b for the cross-section con- 
tributions (2^2 process, 2 — ► 3 process and double-counting term) and the resulting 
matched cross-section for tan/3 = 30, m H ± = 250 GeV, hf = ( m t + m H ± )/^- Note that 
the double-counting term contribution (DC) is subtracted from the sum. 



is collinear with the incoming gluon, and the 2^3 process for large where the 
matrix element should give the correct p±^ distribution. For intermediate p±j < ^Fi 
the double-counting term interpolates smoothly between the two approximations 
giving an overall smooth transition between the collinear and hard regions. We also 
note that the transition region is quite wide ranging in this case from 10 GeV to 
» 100 GeV. 

Restricting ourselves to the experimentally interesting region where 6-quarks 
can be tagged, i.e. p± t b > 20 GeV , we see that the shape of the p^b-distribution 
deviates substantially from the one given by the 2 — * 3 process. In fact, close to this 
lower limit the matched distribution is about twice as large as the 2^3 process. 
In other words, if one only uses the 2-^3 process, the high-px,b region will be 
strongly overestimated compared to the low-pj^i, one. Similarly for rjb (fig. |]b), in 
the experimentally interesting region |^| < 2.5, the shape differs substantially both 
from the leading order and the 2 — > 3 processes, again illustrating that one has to 
use matching in order to get a reliable description of the outgoing 6-quark. 

From fig. |]a it is also clear that the method suggested in [20|| for matching the 
LO and 2 — > 3 processes by making a cut in the p± t b, using the 2-^3 process 
for events with p±^ > P±,cut and the LO process for p±^ < P_i_,cut, rescaled such 
that the total cross-section is given by eq. (|2~5|), does not work in general. Such a 
matching procedure gives the correct cross-section behavior for large p±^ where both 
the leading order process and the double- counting term are negligible, but the draw- 
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Figure 5: Transverse momentum distributions of: (a) the 6-quark from the top decay, 
(b) the hadronic r-jet from — > for tan/3 = 30, m H ± = 250 GeV, and (j,f = 

(m t + m //±)/4. The upper panels show the matched result (solid) compared to the 2 — > 2 
process with parton showers (dashed) , whereas the lower panels show the ratio of the two 
(solid) as well the ratio of the 2^2 process with and without parton showers (dashed). 

back is that the normalization is changed for small p±^ relative to large p±^ since 
the double-counting is only subtracted for p± t i < p±, cut . In addition the kinematic 
constraints represented by the integration limits ( |3.2| ) were not taken into account 
leading to an underestimation of the cross-section, and it is also not guaranteed that 
the differential cross-sections are smooth in all kinematic variables. 

In case one chooses not to tag the accompanying 6-quark the differences between 
the differential distributions from the 2 — > 2 process and the matched ones are not as 
significant. As an example we have chosen to look at the distributions of the 6-quark 
originating from the top quark decay and the hadronic r-jet originating from the 
charged Higgs bosons decay, which are of primary interest when studying hadronic 
r-decays of heavy charged Higgs bosons (see e.g. [|H|). 

The effects of matching on these distributions are shown in fig. [5]. As can be 
seen from fig. |5]a the effects on the shape of the transverse momentum distribution of 
the 6-quark from the top quark decay amounts to 10 — 15% compared to the leading 
order process in the region p^^ < 250 GeV. Similar results are obtained for the r-jet 
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from the decay of the charged Higgs boson, as seen in fig. ^d. These results, except 
for the normalization, are not very sensitive to the choice of factorization scale. As a 
way to gauge the importance of these effects fig. |5] also shows the effect of turning off 
the partons showers on these distributions. From the figure it is clear that with the 
scale for the partons showers set to /ip = (m t + m#±)/4, the effects of matching are 
larger than the effects of the parton shower. However, if one instead sets the scale of 
the parton shower to be (m t + m#±), then the parton shower has a somewhat larger 
impact on the shape of the distributions than the matching. 

At the same time we have found that the effects of matching depend on the cuts 
used in defining the signal. If the cuts are chosen carefully as in [Q the differences 
in the shape of the distributions are diminished. 6 As a consequence it is not pos- 
sible to make a general statement about the need to including matching when the 
accompanying 6-quark is not observed. Instead one has to investigate this from case 
to case if one wants to pin down uncertainties of the size illustrated in fig. ||. 



5. Factorization scale dependence 

Until now we have used what may appear to be a rather small factorization scale 
fip = (m t + rriH±)/A compared to the more conventional choice of fip ~ vl ~ 
m t + iriH±. Of course, in an all orders calculation the factorization scale dependence 
drops out, but when the perturbative series is truncated one is left with a residual 
scale dependence which can be quite large especially in leading order calculations as 
we are dealing with here. However, this formal independence on the factorization 
scale when going to all orders does not mean that the choice is arbitrary. On the 
contrary, as we have already seen, the matrix element for the 2 — > 3 process gives a 
clear indication of where the transition between the collinear and hard regions takes 
place. The important point to keep in mind is that the collinear parton density 
integrates over all p± of the parton up to the factorization scale. In other words, 
we can get a good indication of a proper choice of factorization scale for the leading 
order process from the px^-distribution of the matrix element for the 2 — > 3 process 
shown in fig. |3|. 

In order to be able to extract a suitable factorization scale from the p^-distribution 
it is instructive to compare the matrix element with the double-counting term since 
the latter is based on the collinear approximation. Naively one would expect the 
double-counting term to have a flat p^^-distribution when multiplied with m± b/p±,b 
all the way up to the factorization scale and then drop to zero. However, as is clear 
from fig. |] this is not true. The difference is mainly due to the kinematic constraints 

6 The cuts used in this case are in essence: one hadronic r jet with p± yT > 100 GeV and \ij T \ < 2.5, 
pmiss > -^qq Q e y ; anc j a t i eas t one 6-tagged jet with p±^ > 30 GeV and \rjb\ < 2.5. In addition 
there is also anti-tag against additional 6-jets by requiring that there is not more than one &-jet 
with > 50 GeV and \rj b \ < 2. 
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Figure 6: The differential cross-section da/dp± t b multiplied by m\ b /p±.b for the 2 — > 3 
matrix element and the double-counting term for different factorization scales parameter- 
ized by p = 2//p/(mt + m H ±). 



and a non-constant gluon-density. Comparing the two distributions for different fac- 
torization scales, shown in fig. we see that a suitable factorization scale is obtained 
by requiring that the plateau in the distribution for the double-counting term extends 
to the same p± t b as the matrix element but does not overshoot it. In this way we 
ensure that the size of the collinear logarithms in the 6-quark density are not over- 
estimated. It follows that the factorization scale should be fip ~ {m t + m#±)/4. We 
also see that the "standard" choice \ip = m t + m H ± leads to a large overestimation 
of the collinear logs. 

The p± ^-distribution was also used to find an appropriate factorization scale 
in H and @ when studying the next-to-leading order corrections to the gb — > tH~ 
process. In || the factorization scale was chosen such that the integral over transverse 
momentum of the asymptotic form of the differential cross-section, <^/^P±,fc| asympt = 
S P±,b/ m2 i_ bi gi ves the same total cross-section as the 2 — * 3 process, whereas in [^] one 
compared the asymptotic behavior of the transverse momentum distribution with the 
extent of the plateau in the 2 — > 3 process. Based on these arguments they suggest 
that a suitable factorization scale is fip ~ (m t + m H ±)/6 and [ip ~ {m t + m H ±)/5 
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Figure 7: Integrated cross-sections at m H ± = 250 GeV and tan = 30 as a function of 
the factorization scale parametrized by p = 2/Ltp/ (mj + m#± ) . 

respectively which is very similar to what we find. 

The question of finding an appropriate factorization scale is not particular to 
charged Higgs boson production. In neutral Higgs boson production in association 
with bottom quarks one faces a similar situation where the question of whether to 
treat the incoming 6-quarks as partons or not is also important and one can use 



similar arguments for the appropriate factorization scale to be used. In [27|, f28 
the transverse momentum distribution of 6-quarks in the process gg — > bbh was 
compared to the factorized expectation da/dp± t b oc pj^/m^j,. Using the argument 
that these two distributions should approximately agree up to the factorization scale, 
it was found that a factorization scale of order of m^/10 is more proper than the 
usual (If — 0{mh)- Based on the same argument, a similar conclusion was also 
drawn in case of charged Higgs boson pair-production in association with bottom 
quarks [ 29| . In another approach |30], one instead looked at the distribution 



in Mandelstam-t for the processes bg — > hb and bb — > gh scaled with the Higgs 
boson mass, —tda j ' dt{\f—i /m h ) . Using these arguments the authors argue that a 
small factorization scale [if ~ is preferable. Finally, comparing the NNLO 

calculation of bb — ► h with the LO and NLO ones one |32| finds that both the scale 



dependence and higher order corrections are small when (fi R ,fi F ) ~ {f^hi^hl '4). 

Given that the arguments for choosing the factorization scale to be fiF = {mt + 
m//±)/4 are not very precise we have studied to what extent the results of the match- 
ing procedure changes when varying the factorization scale (in most cases by varying 
it by a factor two up or down) . In order not to confuse the picture we have kept the 
renormalization scale in a s and the running 6-quark mass fixed. 

We start out by considering the factorization scale dependence of the matched 
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Figure 8: The different components of the differential cross-section in p± ^ for different 
choices of the factorization scale parametrized by p = ^ m ^^ ± ) =1/4, 1/2, 1, 2. 

total cross-section. Parameterizing the factorization scale using the parameter p: 

m t + m H ± 

fi F = p (5.1) 

we get the result shown in fig. [7[ Comparing with the scale dependence of the different 
components (2 — > 2, 2 — > 3 and double-counting) which are also displayed we see that 
the scale dependence is substantially reduced after matching. For example, varying 
the factorization scale a factor 2 up or down from p = 0.5 the cross-section changes 
with less than 10%. This reduced factorization scale dependence is expected since 
we are including parts of the NLO corrections (more specifically the real corrections 
from gluon splitting into 66-pairs in the initial state) to the leading order process. 
The same reduction of the scale dependence has also been seen in neutral Higgs 



boson production via 6-quark fusion [ 30 1 . In addition we see that at p > 1 the 



double-counting term becomes larger than the 2^3 process. In view of the fact 
that the double-counting term is included in order to cancel the component of the 
2^3 process that is already included in the leading order process, this marks a 
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Figure 9: Matched ^^-differential cross-sections for different choices of the factorization 
scale > P = {mt 2 Z H±) = V4, 1/2, 1. 

breakdown in the procedure, indicating that the factorization scale is chosen too 
large. 

Next we consider the matching procedure itself. As can be seen from fig. || 
showing the ^-distribution of the 6-quark, the matching procedure still fulfills the 
requirements outlined in section ^| for other factorization scales but the relative im- 
portance of the different contributions varies. However, for the case p > 1, the 
double-counting term is overshooting the 2^3 process in the region p± b < fi F , 
again indicating that this choice of factorization scale is too large. 

Finally fig. shows the relative stability of the matched p^-distribution when 
varying the factorization scale up or down with a factor 2 from our preferred value. 
From the figure we see that the residual factorization scale uncertainty is typically 
small and at most ~ 20% (in the region p±^ ~ (n^t + m #±)/5). This should be 
compared with the difference between the matched distribution and the one from 
the 2 — > 3 matrix element which amounts to a factor ~ 2 at px,b ~ 20 GeV. 

6. Conclusions and Outlook 

In this paper we have presented a new method for matching the gb — > H^t and 
gg — > H^tb processes in Monte Carlo event generators such as Pythia and Her- 
wig. By matching the two processes at hand in a proper way we can combine their 
respective virtues. On the one hand, the gb — > H^t (or 2^2) process includes a re- 
summation of potentially large logarithms of the type a s log(/i^/mfc) which arise from 
the collinear region of phase space where the transverse momentum of the accompa- 



- 17- 



nying 6-quark p± t b is small. On the other hand, the gg — > H^tb (or 2 — > 3) processes 
contains the exact kinematics of the accompanying 6-quark which is important in 
the hard region where is large. 

The matching is done by adding the two different approximations and subtracting 
the common part which otherwise would be double-counted. In other words the 
double-counting term is given by the collinear approximation of the gg — > H tb 
process. By viewing the double-counting term as a differential distribution in the 
kinematic variables of the 2^3 process, it can be used to generate events in the 
same way as any other process, including parton showers and hadronization. These 
events are then subtracted from the sum of the events generated from the gb — > H ± t 
and gg — > H^tb processes. 

This matching procedure leads to a smooth transition between the collinear 
region of phase space, where the 2^2 process dominates, and the hard region 
where the 2^3 process dominates. Looking at the ^^-distribution the difference 
between the matched cross-section and the 2 — ^ 3 process can be as large as a factor 
~ 2 in regions of experimental interest (pjl,& > 20 GeV). Likewise, the pseudorapidity 
distribution of the accompanying 6-quark differs significantly in shape from both the 
2 — > 2 and the 2 — > 3 process in the central region. When looking at the distributions 
of the decay products of the t-quark and the charged Higgs boson the differences are 
smaller, typically ~ 10% in experimentally interesting regions of phase space which 
is similar to the effects of parton showers. At the same time, the differences turn 
out to be sensitive to the cuts used to define the signal and by carefully choosing the 



cuts, as in [Q, the effects can be made negligible. 

Not only does the matching procedure give a better description of charged 
Higgs boson production, in addition it also gives a reduced factorization scale de- 
pendence. The sum of the different contributions to the total cross-section has a 
much smaller factorization scale dependence than the individual parts. Looking 
at the p^-distribution from the gg — ► H tb matrix element and comparing to 
the double-counting term we also get strong arguments for choosing a factorization 
scale \ip — [pit + ^h±)/4 which is substantially smaller than the "standard choice" 
fiF ~ Vs ~ (m t + mH±)- Being conservative we estimate the remaining factorization 
scale dependence by varying the factorization scale with a factor two. Doing this we 
find that the total cross-section varies with about ±10% and that the height of the 
Ri_,b-distribution is also quite stable, varying at most about ±20%. 

The method that we have presented in this paper is not restricted to charged 
Higgs boson production. It can also be applied in other processes where one has 
incoming 6-quarks. The simplest example is gb — > H /± t-production, but it can also 
be extended to include processes with two incoming 6-quarks, such as Higgs boson 
production in association with 6-quarks. It would also be interesting to extend the 
method to also include the remaining NLO corrections to get a NLO normalization 
of the total cross-section as in the MC@NLO method. 
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